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Since  mid-1984,  the  Defense  Mapping  Agency  (DMA)  has  been  computing  predictions  of  pole 
position  (x  and  y)  and  time  difference  (UT1-UTC  =  ^UT)  to  support  the  NAVSTAR/GPS  orbit 
determination  process.  These  predictions  are  based  upon  the  "derived  Earth  orientation  para¬ 
meters"  published  in  the  U.  S.  Naval  Observatory’s  Earth  Orientation  Bulletin  (Series  7) 
and  are  obtained  by  fitting  these  data  to  a  set  of  models  specified  in  the  DMA/GPS  Interface 
Control  Document. 


This  paper  presents  a  review  of  the  models  and  processes  by  which  the  Earth  orientation 
predictions  are  computed  at  DMA.  The  results  of  a  study  to  assess  the  accuracy  of  these 
predictions  are  also  presented.  Finally,  several  proposed  methods  for  improving  the  pre¬ 
diction  accuracies  without  altering  the  basic  models  are  discussed  and  evaluated. 
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THE  DMA/ GPS  EARTH  ORIENTATION  PREDICTION  SERYICE 

JOHN  A.  BANGERT 
DEFENSE  MAPPING  AGENCY 
WASHINGTON,  DC  20305-3000 

SUMMARY 

Since  mil-1984,  fee  Defense  Mapping  Agency  (DMA)  has  been  computing  predictions 
of  pole  position  (x  and  y)  and  time  difference  (UT1-UTC  =  AUT)  to  support  tire 
Navstei/GPS  orbit  determination  process.  These  predictions  are  based  upon  the  "derived 
Earth  orientation  parameters’  pnbhsbed  in  the  U.  S.  Naval  Observatory's  Earth  Orientation 
Bulletin  (Series  7)  and  are  obtained  by  fitting  these  data  to  a  set  of  modeb  specified  in  die 
DMA/GPS  Interface  Control  Document 

This  paper  presents  a  reviev  of  tire  modeb  and  processes  by  which  the  Earth 
orientation  predictions  are  computed  at  DMA.  The  results  of  a  study  to  assess  die  accuracy 
of  these  predictions  are  also  presented.  Finally,  several  proposed  methods  for  improving 
the  prediction  accuracies  without  altering  die  basic  modeb  are  discussed  and  evaluated. 

1.0  INTRODUCTION 

The  GPS  Master  Control  Station  (MCS)  has,  as  one  of  its  primary  functions,  die 
requirement  to  provide  the  GPS  satellites  with  ephemens  and  clock  information  that  can  be 
transmitted  to  users  in  the  navigation  message.  Tracking  das  collected  by  five  monitor 
stations  (MSs)  are  sent  to  the  MCS  and  used,  in  combination  vith  reference  trajectories,  to 
compute  estimates  and  make  predictions  of  ephemens  and  clock  states  for  die  satellites  [1 1. 
These  ephemens  and  cbck  computations  require  die  ability  to  transform  between  an 
Earth-fixed  reference  system  and  a  basic  inertial  system.  Completion  of  this 
transformation  requires  estimates  of  the  Earth  orientation  parameters  (EOP). 

The  EOP  consist  of  the  polar  coordinates  (x  and  y)  and  the  time  difference  UT1-UTC. 
The  polar  coordinates  represent  the  position  of  the  true  celestial  pole  vith  respect  to  a 
reference  point  fixed  to  the  crust  of  die  Earth  (die  Conventional  International  Origin,  or 
CIO  pole).  The  x  coordinate  is  pcsitrre  in  the  direction  of  Greenwich  and  the  y  coordinate 
is  positive  towards  the  vest  [2j.  UT1-UTC  (AUT)  represents  the  difference  between  die 
rotational  time  scab,  UT1,  arid  the  uniform  time  scab,  UTC  (Coordinated  Universal 
Time).  It  can  be  applied  to  UTC  to  recover  time  based  on  die  true  rotetional  rate  of  the 
Earth  (UT1).  Prediction  of  the  EOP  must  be  based  on  extrapolation  of  observations. 

The  Defense  Mapping  Agency  (DMA)  is  required  to  supply  predictions  of  the  EOP  to 


the  MCS  based  on  a  model  and  in  &  format  specified  in  an  Interface  Control  Document 
(ICD-GPS-211)  [3],  This  requirement  is  satisfied  through  the  use  of  an  Eaith  orientation 
parameters  prediction  computer  program  (EOPP),  vhich  is  run  veekly  at  DMA.  This 
program  has  been  providing  predictions  of  the  HOP  to  the  MCS  since  mid- 1 984.  Section  2 
of  this  paper  presents  a  revie  v  of  the  data,  models,  and  processing  scenario  used  by  DMA 
to  produce  the  EOP  predictions.  Section  3  gives  a  summary  of  the  accuracy  of  these 
predictions  from  the  time  the  service  started  in  mid- 1984  to  die  end  of  1985.  A  discussion 
and  evaluation  of  methods  vhich  have  the  potential  to  improve  the  prediction  accuracies  are 
provided  in  Section  4.  Finally,  Section  5  presents  conclusions  vhich  have  resulted  from 
this  study. 

2.0  DATA  PROCESSING 

2.1  EARTH  ORIENTATION  DATA 

The  Earth  orientation  data  used  for  prediction  of  the  EOP  are  provided  by  die  U.S. 
Naval  Observatory  (USNO).  These  data  primarily  consist  of  those  values  of  z,  y,  and 
AUT  published  under  the  heading  of  “derived  Earth  orientation  parameters"  in  the  veekly 
Earth  Orientation  Bulletin  (U.S.  Naval  Observatory  Time  Service  Publication-  Series  7). 
These  values  are  based  on  observations  obtained  solely  from  institutions  vithin  die  United 
States  [4]. 

In  order  to  ensure  the  shortest  possible  delay  betveen  observation  and  data 
dissemination,  USNO  supplies  DMA  vith  the  most  recent  EOP  values  prior  o  publication 
of  the  Series  7  Bulletin.  These  values  are  based  on  the  U SNO  "quick  CORE"  solution  and 
are  available  on  die  General  Electric  (GE)  Mark  III  communications  netvoik  through  die 
USNO  Rapid  Data  Service.  The  "quick  CORE"  pole  values  are  determined  primarily  from 
satellite  laser  ranging  (SLR)  to  die  LAG  EOS  satellite.  The  AUT  values  from  "quick 
CORE"  are  a  combination  of  connected-element  interferometer  (CEI)  observations,  SLR, 
and  optical  observations  [5].  In  late  1985,  die  minimum  delay  betveen  observation  and 
availability  vas  approximately  one  veek.  lAf 

After  publication  of  the  Series  7  Bulletin,  the  "quick  CORE"  EOP  values  in  die  DMA  j 
data  base  are  replaced  by  the  published  values.  These  values  result  from  USNO  "regular !  ' 
CORE"  solutions.  The  "regular  CORE"  pole  values  are  based  primarily  on  SLR 
observations  and  Doppler  tracking  of  the  TRANSIT  navigation  satellites.  The  "regular:  . 
CORE"  AUT  values  are  based  on  a  combination  of  CEI,  SLR,  optical,  and  Doppler 
observations  [5].  Additional  information  concerning  die  USNO  "CORE"  solutions  maybe  — - 
found  in  reference  [6].  Dis 
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2.2  MODELS 

2.2.1  POLE  POSITION 

The  fofloving  modeb  are  used  for  the  pole  coordinates  x  and  y  [3]: 

x  (t)  =  A  +  B  (t-y  ♦  2  Cj  sin  [2ft  (H*)  /Pj]  +  2Dj  cos  [  2ft  (t-^)  /Pj  J  (1) 
y  (t)  =  E  +  F  (t-tg)  ♦  2  G^sia  [211(1-1^X0^]  +  2  Hj^cos  [2ti(t-^>  ]  (2) 

vhere:  ^  »  modified  Johan  dale  (MJD)  of  the  first  day  of  data 
P^Qj  =  fixed  Chandler  period*  435.0  days 
?2>  Q2  =  feed  annual  period  =  365.25  days 
),  k  =  summation  indices  =  1,2 

The  coefficients  A,  Cj,  Dj,  E,  G^,  and  (j,  k  =  1, 2)  are  determined  by  least-squares 
adjustment  In  practice,  the  drift  eras  are  omited  from  the  modeb  so  lhat  the  coefficients 
B  and  F  are  set  equal  to  zero. 

2.2.2  UT1-UTC 

Similarly,  die  AUT  model  [3]  is: 

AUT  (t)  =  I  +  J  (t-tfc)  +  2  Km  sin  (ZtKM^fRgJ  +  2  1^  cos  [2ft  (t-t^^  (3) 
vhere:  =  MJD  of  0  January  of  the  current  year 
m=  summation  index  =  1,2, 3,4 

In  practice,  this  model  b  utilized  by  revriting  die  expression  for  AUT: 

AUT  =  UT1-UTC  =  (UT2-UTC)  -  (UT2-UT1)  (4) 

The  tern  UT2-UT1  represents  the  'seasonal  variation'  and  is  given  by  die 
conventional  formula: 

UT2-UT1- 0.022  sin  pt^Vt^g]  -  0.012  cos  [211(1-^3] 

-  0.006  sin  [2ft(Ht)/R4]  +  0.007  cos  [2n<V'^>yR4]  (5) 

vhere:  R3  =  365.25  days 
R4=  182.625  days 

The  Bureau  International  de  1'  Heure  (BIH)  seasonal  variation  coefficients  are  used  in 
ihb  expression  [7]. 

The  tern  UT2-UTC  represents  the  difference  betveen  UTC  and  the  nearly  uniform 
rotational  time  scale,  UT2,  and  b  modeled  by 

UT2-UTC  =  I  +  J  (t-tjj)  (6) 

Hence,  in  equation  (3),  K3  =  -0.022,  K4  =  0.006,  L3  =  0.012,  and  L4  =  -0.007. 
All  remaining  and  are  set  equal  to  zero  and  the  periods  Rj  and  R3  are  not  used. 
The  coefficients  I  and  J  are  determined  by  least-squares  adjustment.  Detaib  of  the 
processing  scenario  for  AUT  vill  be  given  in  Section  2.3. 


2.3  EOP  PREDICTION  SOFTWARE 

Predictions  of  the  EOP  are  made  through  the  use  of  a  FORTRAN  computer  program 
(EOPP).  This  program  was  designed  and  coded  by  Mr.  Gregory  Arnold  of  DMA  in  1984, 
and  revised  by  the  author  in  late  1985.  The  program  first  extracts  the  most  recent  435  days 
of  pole  position  data  and  the  most  recent  32  days  of  AUT  data  from  the  data  base.  The  435 
days  of  pole  coordinates  span  (approximately)  one  Chandler  period.  The  decision  to  use 
32  day3  of  AUT  data  for  prediction  purposes  was  based  on  a  study  performed  at  USNO 
[  4  ].  This  study  indicated  that  the  smallest  root  mean  square  prediction  errors  (compared 
to  BIH  data),  using  a  model  similar  to  the  one  described  in  Section  2.2.2,  vere  achieved 
using  32  days  of  AUT  data. 

The  EOPP  program  then  performs  an  unweighted  batch  least-squares  fit  to  die  x  data, 
using  the  model  described  in  Section  2.2.1.  Based  on  the  results  of  to  fit,  die  data  are 
edited  so  that  any  point  vhcse  post-fit  residual  lies  outside  three  standard  deviations  is 
tagged.  The  batch  least-squares  fit  is  repeated  using  the  untagged  data.  A  similar 
procedure  is  followed  for  die  y  data.  Since  1  January  1986,  a  method  of  "bias  adjustment- 
has  been  used  to  improve  the  accuracy  of  the  pole  position  predictions.  This  technique, 
which  will  be  discussed  further  in  Section  4,  involves  the  post-fit  adjustment  of  the  x  and  y 
model  bias  terms  (A  and  E)  by  the  value  of  the  residual  at  the  last  untagged  data  point 
This  essentially  forces  die  residual  at  the  last  data  point  to  zero. 

The  reduction  of  the  AUT  data  is  slightly  more  complicated.  First,  die  input  AUT 
values  are  adjusted  for  the  seasonal  variation  by  application  of  equation  (5).  The  resulting 
values  represent  UT2-UTC  (see  equation  (4)).  These  values  are  then  fit,  by  batch  least- 
squares,  to  a  straight  line  (see  equation  (6)).  The  values  are  then  edited  in  the  same  manner 
as  die  pole  data,  and  the  fit  repeated  using  the  untagged  values.  The  seasonal  variations  are 
subtracted  from  UT2-UTC  to  obtain  estimates  of  AUT. 

2.3.1  THE  EOP  PREDICTION  RECORDS 

The  EOPP  program  produces,  as  output,  a  set  of  five  caul-image  records  containing  all 
information  necessary  to  perform  die  EOP  predictions.  The  format  of  these  records  is 
specified  in  ICD-GPS-21 1  and  is  given  in  Table  1 .  In  addition  to  die  coefficients  resulting 
from  die  least-squares  adjustments  (A,Cj,C2,  D^I^E,  Gj,G2,Hj,H2,I,  and  J), 
the  periods  (Pj,  P2,  Qj,  Qo,  Rj,  R2,  %  and  R4),  and  epochs  and  are  given. 
Record  5  also  contains  the  offset  between  International  Atomic  11106  (TAI)  and  UTC. 

The  EOP  prediction  records  are  supplied  to  the  MCS  on  a  magnetic  tape  and  are  also 
placed  in  a  file  on  the  GE  Mark  IH  network  for  immediate  availability.  In  die  near  future, 
esfebhshment  of  a  direct  data  link  between  DMA  and  the  MCS  will  enable  DMA  to  pass  the 


EOP  predxocn  records  to  the  MCS  in  the  most  timely  manner  [  3  j. 

3.0  ACCURACY  OP  THE  DMA  EOP  PREDICTIONS 

3.1  PROCEDURE 

Each  time  the  EOPP  program  is  executed,  the  nevly-comptned  set  of  EOP  prediction 
records  is  appended  to  a  permanent  data  base.  This  data  base  vas  examined  carefully,  and 
a  number  of  duplicate  sets  of  prediction  records  vere  removed.  A  total  of  72  sets  of 
prediction  records  remained  after  this  editing  vas  completed.  The  first  of  these  sets  vas 
produced  in  July  1984  and  ihe  last  3et  used  in  this  study  vas  produced  in  December  1985. 
The  values  from  each  set  of  prediction  records  (see  Section  2.3.1)  vere  inserted  into  the 
models  specified  in  Section  2.2  and  40-day  predictions  of  die  EOP  vere  computed.  In 
each  case ,  da  first  day  of  prediction  vas  the  day  immediately  folio  ving  die  last  day  of  data 
used  to  determine  a  particular  set  of  prediction,  records.  The  predictions  vere  computed  in 
5-day  increments. 

These  predictions  vere  directly  compared  vith  the  actual  values  of  die  EOP  determined 
by  U3NO  and  differences  (in  the  sense,  observation  minus  predicted  value)  vere 
computed.  After  all  72  3ets  of  prediction  records  vere  processed,  root  mean  square  (rm3) 
prediction  errore  at  each  5-day  increment  vere  computed. 

3.2  ACCURACY  OP  POLE  POSITION  PREDICTIONS 

Figure  1  presents  a  plot  of  the  ims  error  in  die  pole  position  predictions  as  a  function 
of  die  age  of  the  prediction.  It  can  be  seen  that  the  x  component  of  the  pole  vas  predicted 
vith  an  rms  error  that  grev  nearly  linearly  from  approximately  0.013  arcsec  5  days  after 
the  last  day  of  data,  to  approximately  0.027  arcsec  40  days  after  the  last  day  of  data.  The 
nns  prediction  error  in  die  y  component  of  the  pole  grev  from  approximately  0.009  arcsec 
to  0.016  arcsec  over  die  same  prediction  interval  Predictions  based  on  all  72  sets  of  EOP 
prediction  records  described  in  Section  3.1  entered  into  the  determination  of  these  rms 
errors.  Note  that  a  prediction  error  of  0.03  arcsec  is  approximately  equivalent  to  an  error 
of  1  meter  at  the  surface  of  die  Earth. 

3.3  ACCURACY  OF  THE  AUT  PREDICTIONS 

Figure  2  is  an  analogous  plot  vhich  gives  the  ims  error  in  predicting  AUT  as  a 
function  of  The  age  of  the  prediction  When  die  nns  errors  based  upon  all  72  sets  of  EOP 
prediction  records  vere  plotted,  it  vas  noticed  that  the  prediction  errors  did  not  increase 
vith  time.  A  detailed  examination  of  the  individual  values  indicated  that  predictions  based 
on  the  first  7  sets  of  EOP  prediction  records  exhibited  unusually  large  prediction  errors 
(20-30  arcsec).  These  large  prediction  erron  appear  to  be  due  to  an  incorrect  treatment  of 


the  seasonal  variation  terras  (see  Section  2.2.2)  in  the  first  7  runs  of  the  EOPP  program. 
Thus,  a  more  realistic  assesment  of  the  AUT  prediction  errors  can  be  made  by  eliminating 
the  first  ?  sets  of  EOP  prediction  records  and  recomputing  the  rms  errors  based  on  the 
remaining  65  sets.  These  errors  are  also  plotted  in  Figure  2.  It  can  be  seen  that  the  nn3 
prediction  error  grovs  from  2.4  msec  5  days  after  the  last  day  of  data  to  7 .5  msec  40  days 
after  the  last  day  of  data.  A  prediction  error  of  2.2  msec  is  approximately  equivalent  to  a 
prediction  error  of  1  meter  at  the  surface  of  the  Earth. 


A  short  study  vas  undertaken  to  investigate  several  techniques  vith  potential  to 
improve  the  accuracy  of  the  EOP  predictions.  It  must  be  emphasised  that  any 
improvements  to  the  EOP  processing  scenario  must  be  consistent  vith  the  models  and 
output  specifications  given  in  ICD-GPS-21 1  (see  Sections  2.2  and  2.3.1).  Given  these 
restrictions,  data  weighting  schemes  vere  given  the  most  attention. 

The  procedure  for  evaluating  these  potential  improvements  vas  as  follovs.  First,  25 
te3t  data  3ets  vere  created  through  extractions  from  die  DMA  data  base  (see  Section  2.1). 
The  first  data  set  commenced  on  MJD  4441 3  (23  June  1980).  The  first  day  of  data  for 
each  subsequent  set  vas  incremented  by  30  days,  so  that  the  last  data  set  began  on  MJD 
45133  (13  June  1932).  Each  data  set  consisted  of  435  days  of  EOP  values.  All  435  days 
of  pole  data  and  the  last  32  days  of  AUT  data  vere  utilized  throughout  testing .  Each  of  die 
25  test  data  sets  vas  processed  by  the  “production"  version  of  the  EOPP  program, 
prediction  records  vere  generated  and  evaluated,  and  rms  prediction  errors  vere  computed 
based  on  the  results  of  all  25  test  runs. 

After  the  addition  of  each  potential  improvement  o  die  EOPP  program,  the  25  lest  runs 
vere  redone.  Again,  the  EOPP  prediction  records  vere  evaluated  and  rms  prediction 
errors  generated.  These  rms  prediction  errors  (at  5-day  intervals)  vere  then  subtracted 
from  the  analogous  prediction  errors  computed  using  the  "production"  results.  The 
differences  betveen  ims  prediction  errors  are  a  measure  of  die  improvement  or  degradation 
of  the  predictions  resulting  from  the  processing  modification.  Each  potential  processing 
improvement  vill  be  discussed  separately. 

4.1.1  BIAS  ADJUSTMENT 

After  die  completion  of  each  EOPP  te3trun,  die  post-fit  residuals  vere  examined.  The 
value  of  the  residual  at  the  last  untagged  data  point  vas  added  id  the  bias  term  (A,  E,  or  I) 
of  die  particular  model  (see  Section  2.2).  This  procedure  forces  the  "nev"  residual  of  the 


last  data  point  to  be  zero.  Thus,  predictions  based  on  the  adjusted  coefficients  should 
begin  unbiased.  Based  upon  the  favorable  results  of  a  previous  study  {8],  this  method  has 
been  routinely  applied  to  the  pole  position  prediction  results  since  1  January  1986  (see 
Section  2.3). 

4.1.2  CONSTRAINED  SOLUTION 

A  very  large  veight  (ic£)  vas  assigned  to  the  last  data  point  This  essentially 
constrained  the  least-squares  solution,  resulting  in  a  very  small  residual  at  the  last  day  of 
data.  Again,  predictions  based  on  this  solution  should  begin  unbiased. 

4.1.3  TIME-WEIGHTED  SOLUTIONS 

Three  different  veighting  schemes,  each  a  function  of  time,  vere  investigated.  The 
general  approach  involved  applying  a  particular  •veighting  function  to  the  la3t  30  days  of 
pole  data  and  (he  last  14  days  of  AUT  data.  These  numbers  vere  chosen  rather  arbitrarily, 
based  on  an  examination  of  production  ran  post-fit  residuals. 

4. 1.3.1  RAMP  WEIGHTING  FUNCTION 

Weights  for  the  affected  EOP  values  vere  determined  by  a  linear  function  of  time.  The 
veights  vere  scaled  30  that  the  last  data  point  received  a  veight  of  1 0  (largest  veight). 

4.1. 3.2  EXPONENTIAL  WEIGHTING  FUNCTION 

Weights  for  die  affected  EOP  values  vere  computed  using  an  exponential  function  of 
time.  Again,  veights  vere  scaled  so  drat  the  last  data  point  received  die  largest  veight  (a 
value  of  10). 


4. 1.3. 3  STEP  WEIGHTING  FUNCTION 

Weights  for  the  affected  EOP  values  vere  set  equal  to  die  same  value  (10). 

4.2  RESULTS 

4.2.1  X-COORDINATE  OP  THE  POLE 

Figure  3  presents  die  differences  betveen  the  rms  prediction  errors  obtained  using  die 
various  processing  modifications  and  the  nns  prediction  errors  obtained  from  die 
production  results.  Each  processing  modification  yielded  a  overall  improvement  in  the 
accuracy  of  the  x -coordinate  predictions.  The  bias  adjustment  and  constrained  solution 
methods  yielded  the  largest  improvements  in  prediction  accuracy.  Use  of  the  exponential, 
step,  and  ramp  veighting  functions  produced  very  similar  improvements  vhich  vere  less 
than  the  improvements  obtained  through  the  use  of  the  first  tvo  methods. 

4.2.2  Y-CQQRDINATE  OF  THE  POLE 

Figure  4  presents  a  similar  plot  of  the  rms  prediction  error  differences  for  the 
y-coordinate  of  the  pole.  Again,  each  processing  modification  resulted  in  an  overall 
improvement  in  prediction  accuracy.  This  improvement,  bovever,  vas  much  smaller  dan 


the  improvement  in  prediction  accuracy  evident  in  the  x-coordinate  results.  In  general,  the 
improvement  in  prediction  accuracy  was  quite  similar  for  all  of  the  tested  techniques. 

4.2.3  AUT 

Figure  5  provides  the  rms  prediction  error  differences  for  AUT.  It  is  evident  that  only 
die  bias  adjustment  method  resulted  in  an  overall  improvement  in  prediction  accuracy. 
This  improvement  is  quite  small  (at  the  half-msec  level).  Use  of  the  remaining  techniques 
resulted,  generally,  in  no  overall  improvement  in  the  AUT  prediction  accuracy. 

5.0  SUMMARY  OF  RESULTS  AHD  COHCLUSIOHS 

DMA  has  been  supplying  sets  of  EOP  prediction  records  to  the  MCS  on  a  weekly 
basis  since  mid- 1984.  These  EOP  prediction  records  are  based  on  unweighted  batch 
least-squares  adjustments  to  the  EOP  data  (x,  y,  and  AUT)  supplied  by  USHO.  The 
accuracy  of  die  EOP  predictions  can  be  assessed  by  comparing  the  predictions  vith  the 
actual  U SNO  EOP  values  available  at  a  later  time.  This  has  been  done  for  72  set3  of  EOP 
prediction  records  produced  from  July  1984  to  December  1985.  The  rms  prediction  errors 
for  the  x  coordinate  of  the  pole  did  not  exceed  0.03  arcsec  (approximately  one  meter  on  the 
surface  of  the  Earth)  over  prediction  intervals  of  40  days.  The  nns  prediction  errors  for  the 
y  coordinate  of  the  pole  did  not  exceed  0.02  arcsec  over  tire  same  prediction  interval.  A 
processing  error  -was  discovered  -which  adversely  affected  AUT  predictions  from  the  first  7 
sets  of  EOP  prediction  records.  When  these  7  sec  vere  omitted  from  analysis,  the  nns 
prediction  errors  for  AUT  were  found  to  be  les3  than  8  msec  over  die  40  day  prediction 
interval 

Several  techniques  with  potential  for  improving  the  accuracy  of  the  ECP  predictions 
vere  evaluated.  These  techniques  included  post-fit  adjustment  of  the  model  bias  terms, 
constraining  die  least-squares  solution  duo  ugh  the  last  data  point,  and  three  weighting 
functions  (a  ramp,  an  exponential,  and  a  step  function)  applied  to  the  kn  30  days  of  pole 
values  and  die  last  14  days  of  AUT  data.  Generally,  all  of  these  techniques  improved  die 
pole  position  prediction  accuracies  during  the  test  period.  The  improvement  ranged  from 
approximately  0.010  arcsec  to  approximately  0.020  arcsec  for  the  x  component,  depending 
upon  technique  and  the  prediction  age.  The  improvement  was  les3  for  the  y  component, 
ranging  from  approximately  0.002  arcsec  to  approximately  0.005  arcsec,  8gain  depending 
upon  method  and  prediction  age.  However,  only  the  bias  adjustment  technique 
improved  the  AUT  predictions.  This  improvement  was  minimal  The  remaining  methods 
generally  produced  no  improvement  in  prediction  accuracy. 

It  is  important  to  note  that  no  attempt  was  made  in  this  study  to  determine  an  optimal 


technique  or  procedure  for  improving  the  EOP  predictions.  Certainly,  other  veighting 
functions  and  vtight  values  could  be  tested.  Experiments  could  be  performed  vhich  vary 
the  amount  of  data  weighted.  However,  die  results  presented  in  this  paper  indicate  that 
techniques  vhich  favor  the  last  days  of  data  lead  to  improved  predictions  of  the  pole 
coordinates  vhen  the  models  specified  in  Section  2.2.1  are  utilized.  These  results  also 
indicate  that  die  same  techniques  applied  to  the  AUT  predictions  are  not  likely  to  lead  to 
improved  predictions  vhen  the  model  specified  in  Section  2.2.2  is  employed. 
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TABLE  1 :  FORMAT  OF  THF  SOP  PREDICTION  RECORDS 
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1 
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61 
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71 
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u re  1:  Prediction  Accuracy  —  Pole  Position 
(72  Prediction  Record  Sets) 
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5  4:  Improvement  in  y  Predictions 
pared  with  "production"  solutions) 

constrained  exponential  step 
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